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An algorithm has been developped, which makes it possible to automatically extract 
trajectories of a large number of particles from fast imaging data, allowing a statisti- 
cal analysis of particles trajectories under various plasma conditions, a better under- 
standing of their influence on plasma properties, and a better characterization of the 
plasma itself. In this contribution, we focus on results obtained in a radiofrequency 
parallel plate reactor, where a large amount of micron-sized carbon dust is produced 
in situ. The use of the rescaled range analysis (R/S analysis) applied to dust particles 
displacements allows decomposing dust dynamic on different time scales. It is shown 
that dust displacement is dominated by collisions on short time scales whereas long 
term behaviour is strongly influenced by large scale plasma fluctuations. 
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I. INTRODUCTION 

During the last decade dusty plasmas have attracted much attention because of their var- 
ious applications, notably in plasma processing technologies. The mechanisms of generation 
and transport of such particles are still extensively investigated. It was shown that, either 
for individual dust grains or for large concentration of dust particles, dynamic and spatial 
distribution strongly depend on a variety of forces acting on thenA Depending on operating 
conditions, equilibrium between these forces varies and is responsible of dust cloud shape 
and dynamic. The investigation of dust particles dynamic can allow obtaining some plasma 
parameters or dust particle characteristics such as their size^. In many cases, oscillations 
are observed which result from a perturbation in the balance between forces acting on dust 
grains. An oscillator equation can usually be derived, showing that the oscillation amplitude 
and period are mainly function of particle mass, charge, and electric field^. Then, working 
with calibrated dust grains, the latter can be used as probes to infer some plasma parameters 
such as the electric field^. On the contrary, the knowledge of plasma parameters makes it 
possible to estimate the mass and/or charge of dust grains by studying their trajectories. In 
both cases, the knowledge of dust particle response to a variation of the control parameter 
is needed. Such studies can be really time consuming and difficult to carry out, especially 
when particles are numerous and/or with complex trajectories. In all cases, it is necessary 
to have access either to a statistical analysis of a large number of dust particles trajectories, 
or to have some accurate dust trajectories on a long duration giving an estimation of the 
overall motion. To perform such an analysis, an algorithm has been developed, which allows 
obtaining dust particles trajectories automatically, from camera measurements showing even 
a large concentration of particles (several thousands) and with a rather simple experimental 
set-up (e.g., no external light source is needed). In this contribution, this powerful diag- 
nostic is used to evidence different types of dust response to sheath oscillations. Paper is 
organized as follow: the description of the experimental setup, ultra fast camera diagnostic 
and operating parameters are given in sec. 2. A brief description of main forces acting on 
dust particles in our plasma configuration is given in sec. 3. Main observations concerning 
the effects of gas injection on dust displacements and collective dynamic are presented in 
sec. 4. Two very different kinds of behavior are underlined in absence and presence of gas 
inlet, which can be explained by the use of the R/S analysis. 
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II. EXPERIMENTAL SETUP 



For investigation of dust transport, macroscopic (tens of /im and larger) carbon dust 
particles were generated in capacitive coupled RF parallel plate reactor working at frequency 
13.56 MHz, with upper electrode being biased and lower one grounded. Vacuum in reactor 
chamber of this device is created by means of rotary vane pump. Plasma was initiated 
in mix of argon and acetylene (C2H2). The motion of dust particles was recorded by an 
ultra fast camera. Besides direct detection and tracking of dust particles, camera observation 
allowed obtaining information about sheath dynamic, such as direct measurements of sheath 
luminescence intensity and thickness. Gases are continuously pumped out, simultaneously 
with gas injection. Mass flow meters control gas injection, which is set independently for each 
component of desired gas mixture. For experiments performed in this study it is necessary 
for dust particles to move approximately in single vertical plane. For this reason, electrodes 
used were rectangular, with length 7 cm, width 2 cm and thickness 0.7 cm. The key element 
of experimental setup is fast camera FASTCAM SAl.l by Photron Ltd. Image recording is 
provided by a 12 bits monochrome CMOS sensor, with a spatial resolution one mega pixel 
up to 5400 fps. To improve resolution, external objective Nikon Micro-NIKKOR 105 mm 
was mounted. With this objective, a resolution of 20 /im per pixel was achieved. Particles 
size depending on operating parameters, these ones must be well chosen to achieve satisfying 
conditions for our experiments. For this purpose, RF power ranges from 10 W to 20 W, and 
pressure of gas mixture (14 % of C2H2 and 86 % of Ar) is 80 mTorr. This range of values 
allows obtaining dust particles with sufficient size to be observable by fast imaging camera. 
For a frame rate of 1000 fps, particle mean displacement between two frames is usually less 
than 1 pixel, which is well suited for tracking. A new diagnostic, still under development, 
based on detection and tracking of dust particles, is used to obtain dust trajectories. It 
can operate in various experimental conditions (low temperature or fusion plasmas) and 
studies (statistical analysis or accurate individual trajectories). A complete description of 
the procedure is discussed in references^. 
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III. FORCES ACTING ON SPATIAL DISTRIBUTION AND TRANSPORT 
OF DUST PARTICLES 

Even in the simplest case (in the absence of magnetic field) of laboratory plasmas, such 
as the following performed experiments, dust particles in plasma volume are subject to the 
action of various forces which confine them in the plasma or drag them outside (see Fig. 2). 
Assuming negatively charged spherical particles, according to the literature^, the charge for 
an isolated dust grain is given by: 

Q d = 4ire r d V d (1) 

with e the vacuum permeability, r d the radius of the dust particle and V d its electrostatic 
potential. All the forces in competition depend on the size of dust particles (15 /xm or more 
in our case), either directly or through the electric charge and then can be distinguished in 
two classes: those independent of the charge (gravitational, neutral drag and thermophoretic 
forces, respectively F g ,F n ,F th in Fig.2 and those directly determined by the charge (elec- 
trostatic force F e i, ion drag force F ion and interaction between particles, which cannot be 
neglected in our case because of the high density of dust)^. In typical capacitive coupled RF 
discharges, like the one used in the following experiments, the most essential forces acting 
on micron-sized dust particles are ion drag force (which tends to confine them), electrostatic 
force and gravitational force (which tend to eject them from the plasma). Note that because 
of the low RF power (typically 10 W), temperature gradient in plasma chamber is small, 
and thus thermophoretic force can be neglected. The following section is dedicated to the 
study of dust particles motion under influence of sheath oscillations. 

IV. INFLUENCE OF SHEATH OSCILLATIONS ON DUST TRANSPORT. 

As already mentioned, normal operation of RF generator used in experiments implies 
continuous inlet and outlet of gas. RF generator can be operated in different regimes when 
gas pumping in and pumping out are switched off so that the amount of gas and pressure 
remain constant. Dust trajectories in presence and in absence of gas inlet exhibit strong 
differences, as illustrated in Fig.3. When there is no gas inlet, particles generally grow at 
a given place and then fall down due to gravity. In the opposite case, trajectories of lots 
of particles exhibit clear tendency for particles to move in the horizontal direction. This 
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direction corresponds to the propagation direction of sheath waves, evidenced by camera 
observations as can be seen in Fig.4. As explained in reference^, these oscillations are directly 
linked to potential fluctuations. By disturbing the force balance, these ones can clearly 
influence dust confinement and transport. This property was used to validate the reliability 
of the tracking algorithm by forcing dust oscillations^. However, in this particular case, low 
frequency oscillations (~ 10 Hz) were mandatory to induce a collective response of dust to 
a forced potential perturbation. In general cases, the frequencies of spontaneous oscillations 
are much higher, of the order of 100-1000 Hz. The fact that sheath oscillations are related 
to gas motion is proved by the observation, that the latter strongly influences the amplitude 
and frequency of sheath fluctuations. To make a comparison, frequency and amplitude 
parameters for two cases only differentiated by gaz injection are adduced in Table 1, with 
light intensity standard deviation being the parameter, characterizing fluctuation amplitude. 
Notably, the amplitude of sheath fluctuations differs in different regions of the sheath namely 
below, above and from the side of the cathode when gas is being pumped in, while when 
there is no gas inlet, deviation is much smaller, indicating that sheath fluctuations are weak, 
and have virtually the same magnitude all around the cathode. In order to describe the 
transport of dust under the influence of a large variety of operation conditions, it is necessary 
to characterize the main feature of particles motion with a simple universal parameter. For 
that purpose, we propose to use the Hurst exponent H. The estimation of H using the 
R/S analysis, is used to evidence any long-period dependance in the temporal behaviour of 
a given phenomenon^. Such an analysis requires the knowledge of successive positions of 
single particles for a time scale statistically large compared to the characteristic evolution 
time of the trajectories. At the moment, the current procedure allows obtaining trajectories 
for about 1000 frames at the best, which is still not fully satisfactory. For that reason, only 
a rather preliminar analysis illustrating the benefits of this technique is presented hereafter. 
For a more accurate analysis, longer trajectories requiring further improvements of the 
tracking algorithm are necessary. The value of H is normalized to 1, which corresponds 
to a purely ballistic trajectory. H = 0.5 characterizes a brownian motion, and values 
between and 0.5 indicate antipersistent phenomena. Fig.5 shows the Hurst exponents H x 
and H y respectively calculated in the directions perpendicular and parallel with respect to 
the cathode, for a single particule under the influence of spontaneous fluctuations seen in 
Fig. 4. Whereas high values of H x and H y point out similar behaviours for short time 
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delays (r < 40 frames), evidencing a classical diplacement between two collisions, a clear 
discrepancy appears for long time delays. The motion in the x direction, perpendicular to 
sheath oscillations, is in between a fractional Brownian motion and a ballistic one, whereas 
dynamic in the y direction is clearly antipersistent. This is a signature of the strong influence 
of the fluctuations on dust transport. These observations are confirmed in all cases where 
long enough trajectories are processed. It proves that the transport of dust, even the massive 
ones, is strongly impacted by fast fluctuations. 

V. CONCLUSION 

In order to investigate dust transport in a laboratory plasma, fast camera imaging was 
performed. Dust trajectories were extracted with an automatic procedure and the advan- 
tages of the R/S analysis for obtaining main features of dust trajectories was demonstrated. 
This technique will be used to study dust transport in different discharges conditions in 
order to determine if plasma fluctuations can be used to control dust transport. 
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TABLE I. Characteristics of light intensity fluctuations in different regions of the sheath. 



Below the cathode Above the cathode From the side of the cathode 



Gaz inlet present Frequency (Hz) 


129.3 


107.7 


110.7 


Standard deviation 


2.7 


19.8 


5.4 


Gas inlet absent Frequency (Hz) 


163.3 


167.8 


163.7 


Standard deviation 


1.4 


1.7 


1.9 




FIG. 2. Representation of the different 

FIG. 1. General scheme of experimental setup. 
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FIG. 3. Comparison of particle trajectories a) in presence and b) in absence of gas inlet. 
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FIG. 4. Temporal dependance of 

sheath light fluctuations along ver- FIG. 5. Hurst exponents corresponding to ver- 
tical (top) and horizontal (bottom) tical (H y ) and horizontal (-Ha;) displacements of 
axes, evidencing vertical propaga- a single dust particle for different time scales, 
tion. 
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